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Abstract 



The Density Dependent Relativistic Hadron Field (DDRH) theory, pre- 
viously introduced and applied to isospin nuclei, is extended to hypernuclei 
by including the octet hyperons. Infinite matter Dirac-Brueckner theory for 
octet baryons and the derivation of in-medium DDRH baryon-meson vertices 
is discussed. From the properties of Dirac-Brueckner interactions it is found 
that hyperon and nucleon self-energies and vertices are related by the ratios 
of free space coupling constants. This leads to simple scaling laws for the 
in-medium hyperon and nucleon vertices. The model is applied in relativistic 
DDRH mean-field calculations to single A nuclei. Free space NA T-matrix 
results are used for the scalar vertex. As the only free parameter the hyperon 
vector vertex scaling factor is adjusted to a selected set of hypernuclear data. 
Spectroscopic data of single A hypernuclei over the full mass range are well 
described. The reduced A spin-orbit splitting is reproduced and found to be 
related closely the medium dependence of scalar and vector interactions. 
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I. INTRODUCTION 



Hypernuclei are unique in providing access to the dynamics of the full meson and baryon 
SU(3) flavour octets. Their study is the natural extension of the isospin dynamics in non- 
strange nuclei towards a more general theory of flavour dynamics in a baryonic environment. 
Obviously, from a QCD point of view hypernuclei as also isospin nuclei are deep in the 
non-perturbative low energy-momentum regime. Hence, a description in terms of mesons 
and baryons should be adequate. Single A hypernuclei, produced in {K~,7i~) or {7i^,K~^) 
reactions on a nuclear neutron are the best studied examples. Their properties confirm 
that adding a unit of strangeness to an isospin nucleus indeed produces a system which, 
to a large extent, follows similar rules as isospin nuclei [0,0. Such observations give strong 
evidence that the strangeness content of a hypernucleus is in fact stored in a hyperon. 
Moreover, hypernuclear spectroscopy indicates the exsistence of shell structures compatible 
with independent (quasi-) particle motion in a static mean-field. The effective potential, 
however, is found to be considerably more shallow than for nucleons. A natural explanation 
for the reduction in depth to about 50% of the nucleon value is provided by assuming that 
the mean-field producing a and u meson fields are not coupled to strangeness. Under these 
"ideal mixing" conditions the meson-hyperon coupling should evolve according to the ratio 
of strange to non-strange quarks in a baryon, i.e. a reduction of vertices by at least a factor 
of RCT,a; ~2/3 is expected for A and E hyperons. However, such a naive quark counting model 
is unable to account for the experimentally observed decrease of the spin-orbit splitting in 
A nuclei. 

Modern approaches to hypernuclear structure are using non-relativistic and relativistic 
microscopic descriptions. Relativisitic mean-field (RMF) theories of Walecka-type have 
been applied successfully with empirically adjusted meson-hyperon vertices. SU(3)/- 
symmetric field theories incorparating chirality or accounting for the quark structure of 
hadrons have been formulated and apphed to hypernuclei. Extensions to still unobserved 
multi-strangeness systems {\S\ >2) have been explored predicting a gain of binding energy 
when adding a few units of strangeness to an isospin nucleus [^p!0|. The production of 
strangelets in ultra-relativistic heavy ion collisions as a new form of hadronic matter has been 
postulated [|lT|. Using SU(3) / arguments nucleon-hyperon and hyperon-hyperon interactions 
in free space [^^-|T^ and in a nuclear environment P^JTB| were calculated. 

In this paper, hypernuclei are described in the Density Dependent Relativistic Hadron 
(DDRH) theory which was introduced previously as an effective field theory for isospin 
nuclei |T7|JI8[] . In DDRH theory the medium dependence of nuclear interactions is described 
by meson-nucleon vertices which are functionals of the fermion field operators. Lorentz- 
invariance, thermodynamical consistency and covariance of the field equations are retained. 
Taking the functional dependence of the vertices on density from infinite matter Dirac- 
Brueckner Hartree-Fock (DBHF) calculations a practically parameter free model Lagrangian 
is obtained once a free space nucleon-nucleon (NN) interaction is chosen. 

A particular conceptual difference to other approaches is the DDRH treatment of non- 
linearities in terms of invariants of fermion field operators. Since the baryon fields are treated 
as quantum fields even in the mean-field limit a well defined class of quantum fluctuations 



with non- vanishing ground state expectation value is taken into account |T8| , p!9| . Dynami- 



cally, they contribute to the Dirac equations as rearrangement self-energies describing the 
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static polarization of the medium. In standard relativistic mean-field (RMF) theory non- 
linearities are attributed to higher order self-interactions of meson fields [Q. In mean- field 
approximation, mesons are treated as classical fields and fluctuations around the classical 
field configurations are neglected by definition. In bulk quantities, as for example total 
binding energies, the differences of the DDRH and the RMF approach are hardly detectable 
because the DDRH rearrangement self-energies are cancelled exactly in extensive thermo- 
dynamical quantities [0. But in single particle quantities like separation energies, wave 
functions and density matrices the differences become visible |]18[. DDRH coupling con- 



stants in asymmetric matter and modifications from vacuum polarization have 
been investigated. 

The extension of DDRH theory to strange baryons is discussed in sect. |T|. The theoretical 
formulation is kept general allowing to include the lowest SU(3) baryon and meson octetts. 
By practical reasons, however, only strangeness-neutral meson fields are taken into account 
at present. As the central theoretical result we derive in sect. 11 B scaling laws for in- 
medium hyperon vertices given by an almost density independent renormalization of the 
nucleon vertices through the ratios of free space coupling constants. In phenomenological 
RMF approaches P,^p3|, a similar scaling ansatz for the meson-A vertices is used but here it 
is obtained theoretically. Equations of motion and the Hartree mean-field limit are derived. 
In sect. |TTl|a reduced model appropriate for single A nuclei, is introduced. Since strangeness 
carrying mean-fields can be neglected (they are of the order 0{1/A), A the mass number) 
the mean-field equations are considerably simplified. In the applications the a coupling is 
taken from a theoretical NA T-matrix [ll2| , p^ while the u coupling is determined empirically. 
DDRH mean-field results for hypernuclei are presented in sect. and compared to data and 
conventional RMF calculations. Spectroscopic data are well described, reproducing also the 
reduced spin-orbit splitting in A nuclei. The paper closes with a summary and conclusions 
in sect. 0. 



II. DENSITY DEPENDENT HADRON FIELD THEORY WITH HYPERONS 

A. The model Lagrangian 

The derivation of a symmetry-broken physical model from a SU(3)/ Lagrangian has 
been excercised e.g. in ref. [Q. However, in order to describe nuclear structure phenomena 
one finds that most of the SU(3)/ structures are actually not contributing. The reason 
is obvious because parity conservation inhibits the appearance of condensed pseudoscalar 
fields in nuclei. Hence, neither of the 0~ meson fields contributes directly to a hypernuclear 
calculation, except through exchange interactions from antisymmetrization. From the 1~ 
vector meson octett condensed isoscalar uj and isovector p meson fields will evolve. In a 
system with a large fraction of hyperons also condensed octett K* and singlett $ mesons 
fields can appear. However, an apparent shortcoming of a pure SU(3)j approach is the 
missing of scalar mesons and, hence, the absence of a binding mean-field. A satisfactory 
description of the O"*" meson channels, e.g. in terms of dynamical 2-meson correlations |lT2| , p4| , 
is an unsolved question. 



3 



Here, we follow the line of relativistic mean-field theory |^ and restrict the model to the 
degrees of freedom which are relevant for the nuclear structure problem. In practice, we use 
the DDRH Lagrangian which is extended in the baryon sector by including the lowest 
S = — 1 (A, E) and S = — 2 (S) baryons. We introduce the flavour spinor "i/p 

= (^jv,^a,^s,^h)^ (1) 
being composed of the isospin multiplets 

where ipi are Dirac spinors. The full Lagrangian is structured in an isospin symmetric way. 
In the exchange particle sector the isoscalar a, as (= scalar ss condensate), to and meson, 
the isovector p meson and the photon 7 are included. This leads to the Lagrangian 



= i E {d,'^^d^<l>. - ml<l>^) - \ E (IfW^ - mlA^-A (3) 

i=(T,as K=u),<f>,p,'y 
+^Ft,^(^F, ^F)^FCrs - ^(^F, ^Fhp^F(l)^ " e^pQlp^FAf", 

where Cb and Cm are the free baryonic and mesonic Lagrangians, respectively, and the 
diagonal matrix M contains the free-space baryon masses. The meson-baryon interactions 
are contained in Cint- 

= d,Ai^^ - S.AW (4) 

is the field strength tensor of either the vector mesons {k = u, 0, p) or the photon [k = 
7). In eq. (0) contractions of the field strength tensors are abbreviated as F"^ = F^^F^^ 
etc.. Q is the electric charge operator, as and $ meson fields are included mainly for 
reasons of completeness. The corresponding classical, condensed fields will be important 
only in hypermatter with a significant strangness content, i.e. for hyperon and nucleon 
fractions of comparable order. For single A nuclei they should be negligible. Lagrangians 
of a similar structure, but with constant meson-baryon vertices, are used successfully in 
relativistic mean- field calculations of hypernuclei, see e.g. 

An isospin symmetric interaction is obtained with vertices chosen as: 



P ^ _P . V -V -B fori?,i?' = Ar,A,S,S, (5) 

J- p I —i^pBOBB', ^ pB — pBT 
/ BB' 

where are the isospin Pauli- matrices. In DDRH theory the vertices F^s = ^aBi"^, ^) are 
taken as functionals of the baryon field operators [|1^,|TH]. 
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B. Baryon-baryon vertices from Dirac-Brueckner theory 



In order to understand the subtleties of including baryon-baryon (BB) correlations into a 
field theory of a higher fiavour content we briefiy sketch the derivation of the vertices F from 
Dirac-Brueckner theory. The main outcome of the discussions is that nucleon and hyperon 
dynamics should be related to a good approximation by simple scaling laws. 

A Lagrangian of the type as defined above leads to a ladder kernel V^^ given 
in momentum represention by the superposition of one boson exchange (OBE) potentials 
Vf^'(q',q). The latter are Lorentz-invariants 

Vr'{q',q) = 9aB9aB'D^{t){u^-{q')K^u^-^{q)) ■ {u^'2{-q')n'^u\-q)) (6) 

where t=(g' — g)^ is the 4-momentum transfer, Kq, denotes the Dirac and fiavour structure 
of the vertex with bare coupling constants gaB and QqB' for baryons B and B' belonging 
to different isospin multiplets. The Dirac spinors are indicated by u^{q). Contractions 
over Dirac and fiavour indices are indicated by the notation ■ k". Working in the BB' 
center-of-momentum frame q and q' denote the relative 4-momenta in the in- and outgoing 
channels (-B3, B'^) and (-Bi, B2), respectively, where -q=(qorq)- 

Solving the Bethe-Salpeter equation with the kernel V^^ = J2a ^ 

n^^\q',q) = V^^\q\q) + J dk V^^' {q' , k)GQp{BB'){k, q)n^^' {k, q) . (7) 

leads to the in-medium interactions 7?. = {kf{q')k2'{—q')\TZ{kp,kp)\kf{q)kf'{—q)). The 
Pauli-projected intermediate two-particle propagation, denoted by GQp{BB'), introduces an 
intrinsic dependence on the Fermi momenta {kp,kp ). Dependencies on the (conserved) total 
center-of-mass energy s = {k^ + fcf )^ are implicit. Evaluating eq. (|^) with self-consistent 
in-medium spinors, including the self-energies S'^(A;), introduces additional medium depen- 
dencies. 

In structure calculations the G-matrices TZ^^ are required in the nuclear matter rest 
frame rather than in the 2-body cm. system. In practice, the transformation is achieved 
by projection on the standard set of scalar (S), vector (V), tensor (T), axial vector (A) and 



pseudo scalar (P) Lorentz invariants p6|-p8|]. For our purpose, however, a more convenient 
representation is obtained by forming appropriate linear combinations of the so obtained co- 
efficients and to map the set (S,V,T,A,P) onto the vertices of OBE interactions, eq. (^. This 
allows to express the {BB') G-matrices, eq. (|^), in terms of renormalized OBE interactions, 

7^f^'(g',g) = z^'''{s,t,u\k^kf)g^sgaB'D^{t){u''^{q')K^u''^{q)) ■ {u^'^{-q')K'^u^'^{-q)) (8) 

The vertex invariants have been decomposed further into boson propagators Da{t) and 
renormalization coefficients z^^ which both are Lorentz invariants. From this representation 
it is apparent that correlations are shifted into the vertex factors z^^ . In principle, they 
may depend on the full set of Mandelstam variables s, t and u and the Fermi momenta 
of baryons B and B'. However, most of the t dependence is already accounted for by 
the meson propagator Da{t) and considering the mild variation of DB G-matrices on the 
center-of-mass energy the z-coefficients can be expected to depend mainly on the Fermi 
momenta. If necessary, the u dependence can be removed to a large extent by adding 
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a term proportional to Da{u), i.e. introducing antisymmetrization explicitly. Note, that 
antisymmetrization effects contribute only to states within the same isospin multiplet, i.e. 
B = B' . The self-energy of baryon S, however, includes contributions from all multiplets: 

(M^(A;)fi:°M^(A;))Ef = 

E / A (Tr((A:^g^'|7^^^>^g^V|'(g)) + 5^^,(fc^g^'|7^^^'|g^A:^V|'(g)) (9) 

B' 

where the space-like integration extends over the Fermi spheres Kp of baryons B' with in- 
medium (positive energy) propagators Sp (g) and 'kp = {kp, kp) denotes the set of nucleon 
and hyperon Fermi momenta. 

A particularly appealing aspect of eq. (H) is that the z-coefficients can be considered as 
medium- dependent renormalization factors of the OBE vertices. Exploiting the fact that 
the SU(3) / isospin multiplets are not mixed by strong interactions we are allowed to assume 
separability 

zr(fc,g|fc|,fcf')^.f(fc|).f(fc|') (10) 
and to neglect the (weak) residual momentum dependence. As a consequence, 

9aBz!^'''{k,q)g^B' ^T^B{k^)TaB'{k^') , (11) 

where 

raBik^)=gaBS^{k^) (12) 

defines the renormalized in-medium vertices in ladder approximation. We introduce the 
antisymmetrized condensed Dirac Hartree-Fock (DHF) meson fields 

(P^{k\kF,r)=Y,r^B'D^{0) [ ^^dqTr{{u''\q)K^u'''iq))) (13) 

B' 

+ J2jj^B^^ f'^i^^i^BDf,{k - q){u^{q)KaU^iq)) 



where /^^ denotes the Fierz matrix [0,^] and F = (F^, F^). With our choice of momentum 
independent, global vertices, eq. (|]) takes then the approximate form 

^aik\k^, kp) - raB{k^)(pa{k\k^, k^ , F) (14) 

This equation establishes the link to the DDRH Lagrangian: In DHF approximation self- 
energies of the same structure are obtained from eq. (|]). 

In order to derive a self-contained model we apparently have to introduce a "renormal- 
ization" scheme. This is achieved by choosing symmetric hyper matter, i.e. kp = = kp- 
Writing down for that case eq. for nucleons and hyperons explicitly one finds 



^avikp) = TaN{kF)^q^^^\k=kF,k^=ky,- (15) 
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which is exact in Hartree approximation. The above relation is the central result of this 
section. The obvious medium dependencies introduced in eqs. (§) and (Q) by the external 
baryon lines are eliminated such that the intrinsic medium properties of the underlying 
interactions are projected out. This is seen more clearly by considering the diagrammatic 
structure of the DB self-energies p6| , p7| . The leading order (Hartree) contribution is given 
by tadpole diagrams and from a perturbation series expansion in the bare coupling constants 

9aB 

Rl = %^'-^{l + 0{l-^)) + ... (16) 

9aN My 

where the reahstic case Qoy < QaN is considered. Hence, the i?^ are expected to be state- 
independent, universal constants whose values are close to the ratios of the bare coupling 
constants. For asymmetric matter with a hyperon fraction = ^ <^ 1 a corresponding 
diagrammatic analysis shows that asymmetry terms are in fact suppressed because the asym- 
metry correction is of leading second order 0{{j^C,yY)- Thus, even in a finite nucleus where 
C^Y may vary over the nuclear volume, we expect = const, to a very good approximation. 

In fact, these results agree with the conclusions drawn from the analyses of single hy- 
pernuclei in the past. In the present context, eqs. (^) and ( [T6| ) are of particular interest 
because they allow to extend the DDRH approach in a theoretically meaningful way to 
hypernuclei using the results available already from the previous investigations of systems 
without strangeness. In the applications discussed below the nucleon (Hartree) scalar and 
vector vertex functions T^^^N^kp) of |jT8| will be used. The hyperon scaling factors are 
treated as phenomenological constants to be determined empirically. In the scalar channel 
information on available from recent calculations of the Jiilich group for the free NA 
T- matrix p!^j2^ is taken into account leaving essentially the ratio R^ / R^ for Y = A as the 



only free parameter. 

We close this section by remarking that eq. (|T^) actually defines a set of quadratic 
equations for the vertices T^b, as seen immediately when inserting eq. ([T3|) into eq. (|14|). 
Vertices derived in this way would be appropriate for DDRH calculations in Dirac-Hartree 
Fock approximation. From the a — u model it is known that DHF and relativistic Hartree 
calculations give almost undistinguishable results for properly adjusted parameters f^. In 
the following we take advantage of that observation and, as in |T8[ , restrict the calculations 
to the Hartree case only. This corresponds to determine the DDRH vertices by expressing 
the $Q, fields on the right hand side of eq. (p!^) in Hartree approximation. 



C. The equations of motion 

In DDRH theory the above results are embedded into a relativistically covariant and 
thermodynamically consistent field theory from which the vertices are retrieved when evalu- 
ated in mean-field approximation. As discussed in [T^ the kp dependence of the DB vertices 
is expressed in terms of Lorentz-scalar (products of) bilinear forms p of the baryonic field 
operators "^p- This provides a unique mapping of the medium dependence onto frame- 
independent Lorentz scalar quantities. The external Dirac structure of the vertices is fully 
determined by the Lorentz character of the meson field. The intrinsic density dependence 
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must be deduced from microscopic calculations as discussed in the previous section. As an 
obvious generalization of the ansatz used in |jl8| the DDRH vertices are expressed here as 



r = r^ij(p„B(^F,^F)), 
a = (7,as,uj,(f), p B = N,A,T,,E 



(17) 



where paB denotes a Lorentz-scalar combination of the baryon field operators. 

By definition the DB vertices ^aEi^p) are c-number valued functions of the Hartree or 
Hartree-Fock expectation value of p^E- From a general theoretical point of view the DDRH 
vertices TaBiPaB) are not necessarily restricted to this particular sub-class of diagrams. 
Formally, a projection onto DB correlations is defined by the mapping |T^ 

(•OO 

TMPaB)= T^i{p^iW^i-PaB)dp^i . (18) 

Jo 

In the following a straightforward extension of the vector density dependence (VDD) 



prescription of ref. [|l^ will be used. This corresponds to the ansatz 



(19) 



and chosing (^B 



aB 



B'B" 



u^^bn5b'n^b"n with a four- velocity (see 



ra). We thus find 



PaB = yjj^j^^ and by means of eq. (|T8|) a Taylor series expansion of the DDRH vertex in 
terms (of the modulus) of only the respective baryon four-vector current is obtained. Ex- 



actly that choice of i?"^ is a practical implementation of the results obtained in sect. p^lB 



This shows that the DDRH approach in fact corresponds to expressing many-body correla- 
tions by an expansion of vertices into baryon n-point functions chosen such that in ground 
state expectation values correlation diagrams of the fully interacting theory are cancelled by 
compensating terms in the DDRH expansion. 

As pointed out already in ref. [|18| the most important difference of DDRH and RMF 
theory are contributions from rearrangement self-energies to the DDRH baryon field equa- 
tions. Rearrangement self-energies account physically for static polarization effects in the 
nuclear medium, cancelling certain classes of particle-hole diagrams Due to the ad- 

ditional strangeness degree of freedom the structure of these rearrangement self-energies is 
much more complex than in the purely nucleonic DDRH theory. The variational derivative 
of Cint now leads to 



6Cr, 



int 



dC 



int 



+ 



E 



dCint Sp, 



aB 



B=N,A,i:,E 



dpaB S^j 



(20) 



With = ^ one finds 

OpaB 



aB 



- dr^{^F,'^F),,, ^„ ^ dr^sipaB 



B 



dpa 

dTaB^PaBj ^a '-p 



dp. 



aB 



dpa 



B 



KJb^ a = w,p. 



(21) 
(22) 
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where and pf are the vector current operator and the scalar density operator of baryon 
type B, respectively. The rearrangement self-energy thus is given by 

f.{r) ^ f.{r)^^, ^23) 



B 



(24) 



where the sums in eq. (0) are those appearing in eq. (pOl). The usual self-energies are 
given through 



(25) 
(26) 



where the f$ are those defined in eq. (^. Thus, the total baryon self-energies are finally 
obtained as 



(27) 



Here, the Fq, are diagonal matrices containing the flavour dependent vertices. However, in 
structure the baryon field equations remain unchanged 



= 



(28) 



D. Mean-field theory 

A solvable model is obtained in the Hartree mean-field approximation which amounts to 
assume that products of fermion operators are normal ordered with respect to the Hartree 
ground state |0), given by a single slater determinant of occupied fermion states. Expecta- 
tions values with respect to the Hartree ground state will be abbreviated as = (0| O |0). 
In the Hartree approach the vertex functionals TasipaB) can be treated in a particularly 
simple way. Applying Wick's theorem one gets [jl8| 

{^aBiPaB)) = TasiPaB), PaB = {Pob) (29) 

which brings the originally highly nonlinear field equations into a tractable form. Corre- 
spondingly the rearrangement contributions are obtained as 

dVaBiPaB) \ _ dTgBiPaB) .^qx 
dpaB I dpaB 

In the approximation as static classical fields the meson field equations reduce to 
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(-V2 + my$„= ^Mp^.B)pf, a = a,as (31) 

B=Ar,A,E,H 

(-V2 + my$;^= Y: ^MP'^.b)^', a = u,c^ (32) 

B=Ar,A,E,H 

(-V2 + mj)p-^= ^ f,B(PpiJ)J^^ (33) 
b=n,a,j:,e 

the Dirac equation for the baryons remains the only equation of motion for an operator field, 
with dynamics given now by static but density dependent self-energies 

- S^(p)) - {m - S,) 
In finite nuclei, where p = p(r), S(p) = S(r) depends on the spatial coordinates. 



0. 



(34) 



III. THE A iV MODEL 

In order to test the scaling relation derived in sect. |l| on existing hypernuclear data a 
model with A-hyperons and nucleons interacting only by non-strange mesons will be dis- 
cussed. Naive quark counting suggests that even for A-A interactions the strange mesons 
only contribute about 10% of the interaction strength. Keeping in mind the large uncer- 
tainties in the hyperon-nucleon and hyperon-hyperon interactions the exchange of strange 
mesons can be safely absorbed in the A-cr and the A-u vertices. We will use the extended 
VDD prescription introduced in sect. II C with a density independent N-p couplings, which 



leads to rather satisfactory results for isospin nuclei 



A. A meson interaction 

Due to the simple interaction structure of A hyperons - the As are isoscalar and elec- 
trically neutral and thus couple neither to the p meson nor to the Coulomb field - the 
investigation of the A-nucleon interaction becomes rather transparent in this model. As- 
suming that the a^- and the 0-meson are pure ss-states and therefore mainly couple to 
the strange quark (due to OZI supression), they will have no significant effect in single-A 
hypernuclei and can thus be safely neglected. 



According to sect. II B the density dependent vertices for the A-N model are given by: 



^uA = Rlu ■ '^LuNiPA) (35) 



Pa = Pa A = puA = yjf^j^^, 



where R„„, from now on denotes R„, ,. The values for R„,„ will be determined in sect. IIII C 



The parameterizations of the nucleon-meson vertices Tij^{p]\f) are taken from Numeri- 
cally a fit with a second order polynomial in fc^;' to the vertices derived in |^T| from nuclear 
matter DBHF self-energies is used. 



10 



B. A rearrangement dynamics 



The considerations of section |IIl7\| define the dynamics of the A-nucleon system, i.e. the 
usual and the rearrangement self-energies can now be specified. Going to the nuclear rest 
frame, the arguments of the vertex functional, p^B, are now defined as in eq. (|19D with 



B 



5r 



p 









(36) 



and it^ = (1, 0, 0, 0). This leads to the rearrangement self-energy 



B=N,A 



or,. 



-a 



_dpaB dp^B 

where the nucleon and A parts are given explicitly by 



(37) 



u 



B 



(38) 



and B = N, A vertices depend intrinsically only on their own densities, as derived in 
sect. ITTBl 



C. The vertex scaling factors 

A consistent extension of the DDRH theory to strangeness requires to use vertex func- 
tionals from DB self-energies calculated in the complete octet sector. However, since such a 
full scale calculation is neither available nor feasible under the present conditions, we choose 
a semi-empirical approach combining existing theoretical information on the A — a vertex 
with a phenomenological description of the u) vertex scaling factor. 



Actually, an extension of the Bonn potential to the free A^-A system already exists ||12| , |24 
but DB calculations are pending. We use the extended Bonn A potential as guideline to 
determine the relative couplings and R^. This is consistent with the approach used in the 
isospin sector because the nucleonic DDRH parameters also were derived from the Bonn-A 
potential [|l^. Because the DB interactions include highly non-linear and non-perturbative 



correlation effects the quark model reduction factor R^ = 2/3 is not expected to be adequate 
for hypernuclear structure studies. 

Clearly, the final decision on the permissible {Ra, Ruj) pairs is obtained from a compari- 
sion to data. In fig. |l] the deviations of calculated and measured A single particle spectra 
are shown. Varying freely the R^ and R^j scaling factors, DDRH single particle energies 
for 208pb^, 89Ya, 51Va, 4°CaA, ^^SIa, ^^Qa, ^^Ca and ^Bca were compared to data deduced 



from (tt, K) [p^-|36| experiments. With precise measurements resolving spin-orbit doublets 
to high accuracy such a procedure would, in fact, allow to fix both the scalar and vector 
scaling factors unambigously because the centroid and splitting energies are determined by 
difference and sum of the scalar and vector mean-field components, respectively. 

Unfortunately, under the present experimental conditions doublets are not resolved ener- 
getically. Typically, spin-orbit splittings are deduced rather indirectly, e.g. by a line shape 
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analysis ||3^ which only allows to set constraints on the energy splitting of spin-orbit part- 



The probably best known case is the (IP3/2IP1/2) splitting in ^^C\. A single data point 



ners. The consequences of these experimental uncertainities for a theoretical analysis are 
clearly seen in fig. |1|: The distribution is characterized by a sharp deep valley extending 
between R^) = (0.1, 0.1) and {R^, R^) = (0.8, 0.9) and any (i?o-, -Rtj)-pair in this region 
would describe the data almost equally well. The quark model value pair R„ = R^ = Rq 
{Rq = |) is seen to be located at the ridge of valley. These results clearly illustrate the 
necessity of high resolution measurements which may become possible in the near future 
with a new generation of detection systems. 

53/2-LPl/2j 

is available from measurements at Brookhaven in 1981 ||3^ from which AE = (0.36 ± 0.3) 
MeV was deduced. In the x^'Pi'ocedure this constraint only excludes extreme values of 
Ro- and R,^ that are already ruled out by the systematics of hyperon binding-energies (see 
also sec. |1V C[) , anyway. Despite the large uncertainties the data seem to favour a small A 
spin-orbit splitting. Very recent measurements at AGS/E929 apparently confirm this 
conclusion. 

In order to remain as close as possible to the microscope DDRH picture we use 

R„ = 0.490 (39) 



which was extracted by Haidenbauer et al. for a sharp a meson of mass mo-=550 MeV. 



Because in |T^J2^ the scalar meson channels were described by the correlated exchange of 
pion and kaon pairs the scalar coupling also includes admixtures of a 0"^ ~ ss field being 
relevant for the A couplings. Since theoretical values for the u vertex are not available R^ is 
treated as a phenomenological parameter. From Fig. |^ and the above value of R^r one finds 
immediately 

R^ = 0.553 . (40) 



In a constant coupling RMF model by Ma et al. [^] for the same value of a rela- 
tive u coupling R^=0.512 was obtained. Considering the quite different DDRH interaction 
structure the deviation is only apparent and, in fact, a surprisingly good agreement can be 
stated. Moreover, our relative couplings are also consistent with bounds on hyperon-nucleon 
couplings extracted from neutron star models P,^^. 

Very likely, most of the deviations of R^- and R^^ to the quark model value of 2/3 are 
caused by the highly nonlinear contributions from the dynamically generated a and ag 
exchange channel in [|12|,^. These genuine many-body effects superimpose additional con- 
tributions from the explicit SU(3)/ symmetry breaking and uo — cf) octet-singlet mixing on 
the fundamental strong interaction level. Obviously, all these effects cooperate in the same 
direction, namely to produce deviations from the limiting values predicted by exact SU(3) 
symmetry. Apparently, the attempt to represent the rather involved dynamics of the scalar 
channel by a single meson of sharp mass implies an effective vector field of compensating 
repulsive strength. As a consequence, neither of them can be expected to resemble the 
properties of the respective bare physical meson. 

In SU(3)-symmetric models explicit symmetry breaking must be introduced in order to 
reproduce hypernuclear spectra, e.g. by a symmetry breaking term in the Lagrangian of 
in the generalized chiral SU(3) cr-model [H] or by means of vertex scaling factors as in the 
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quark- meson-coupling (QMC) model 0. Because of the close similarity of the QMC and 
our approach it is instructive to compare to the results of Ref. [0: For R„ = Rg, has to 
be rescaled by a factor of 0.93, while R^^ = Rg requires to multiply R^j by 1.10. As can be 
extracted easily from fig. [l|, for R^^ = Rg we would get R^^ = Rg ■ 0.97, and, vice versa, for 
Rw = Rq we find R^r = Rg ■ 1.15. A speciality of the QMC is a tensor-coupling term arising 
from the quark-structure of the baryons which keeps the A spin-orbit splitting extremely 
small although the SU(3) couplings are rather large. 



IV. RELATIVISTIC HARTREE DESCRIPTION OF SINGLE A HYPERNUCLEI 



Relativistic DDRH Hartree theory and applications to isospin nuclei were discussed in 
great detail in ref. [0] and the references therein. Here, we present DDRH results only 
for single A hypernuclei. The numerical realization follows closely ref. 0] namely the 
meson fields are described by eqs. (^T])-(^3|) and baryonic wave functions are obtained from 
The nucleon-meson coupling functionals are those of The model parameters 



eq 



are compiled in table |. 



A. Density dependent A vertices in finite nuclei 

Numerically, DDRH calculations rely on baryon-meson vertices taken from infinite matter 
DB interactions which are applied to finite nuclei in local density approximation (LDA). 
The success of DDRH theory in describing isospin nuclei is closely related to the saturation 
properties of nuclear densities from which it is clear that inifinite matter conditions are 
approached gradually with increasing mass number. Under such conditions the LDA is 
likely to be a rather reliable approach. For light nuclei or, as in single A nuclei, small 
fractions of a specific baryon species with respect to the bulk components effects from the 
finite size and finite particle number could limit the applicability. 

The variation in the effective coupling strengths over the mass table is illustrated in the 
upper graph of fig. ^for the single A hypernuclei ^Bca, ^^Oa and ^''^PbA. The lower graph of 
fig. displays the vector densities for IS1/2A states. In the light nuclei the coupling decreases 
rapidly towards the nuclear center while in lead the DDRH vertices are almost constant. The 
behaviour follows closely the density distributions of the Is A states. Their radial extensions 
are determined by the size of the mean-field produced by the nuclear core. Of particular 
interest is the variation of the vertex functionals with nuclear mass. Apparently, the global 
density dependence of the infinite matter DB couplings transforms in finite nuclei effectively 
into a pronounced mass number dependence of the DDRH vertices. From these results it 
is clear that a complete test of the medium dependence is only obtained in calculations 
over a wide mass range. Light hypernuclei will be most appropriate to study those vertex 
parts depending explicitly on the density while heavy nuclei mainly provide information on 
hyperon interactions at a saturated density. This also points to possible limitations of the 
present model: the use of LDA vertices may lead to uncertainties in light nuclei where a 
particular sensitivity on the transition to free space conditions appears. 
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B. Structure of the A mean-field 



An appropriate way to understand hypernuclear dynamics and to compare to other 
calculations is to consider the Schroedinger-equivalent potentials given in lowest order by 
the difference and (the gradient of) the sum of the relativistic scalar and vector mean- fields 
for the central and the spin-orbit potentials, respectively . 

It is obvious, that the scalar and vector fields scale according to the A-meson coupling. 
Similar to other approaches and in agreement with empirical anaylses the Schroedinger- 
type DDRH A-nucleus potential is reduced by a factor of 0.35 to 0.4 compared to the 
nucleon potential. Results are displayed in tab. || and fig. |^. Since binding energies are 
reduced accordingly already the wave functions of the deepest bound A states are spread 
over a large part of the nuclear volume resulting in a sensitivity to the complete surrounding 
density structure. Tab. showing the rms-radii of A, neutron and proton states in the 
hypernuclei ^°CaA and ^°^PbA, illustrate this effect. In fig. ^ the conventional central and the 
rearrangement central A potential for light to heavy nuclei are displayed. It is clearly seen 
that the rearrangement polarization effects are most important in the surface dominated light 
nuclei. Still, the rearrangement self-energies play only a minor role for single A hypernuclei 
since they are weighted by the A vector density (see eqs. (|21|), (p2D) which is obviously fairly 
small. The more important effect of the density dependent treatment arises in hypernuclei 
through the nuclear core creating the A's mean-field potentials. The density dependence of 
the nucleon vertices modifies the core density distribution over the whole nuclear volume 
thereby directly affecting the A mean-field. Since the potential shape and strength 



is reflected in single particle energies and wave functions a A acts as an external probe 
providing a global measure of the core properties. 



C. Single particle states 

Hyperon single particle spectra for IS"! = 1 hypernuclei can be seen as a very clean 
fingerprint of this nucleus, since, as discussed in the last section, they are almost undisturbed 
by many-body effects. Besides the bulk structure, which contains information on the mean- 
field, i.e. the nucleonic density distribution, the spectra yield information also on other 
correlations of the baryonic interaction, carried by the fine structure. 

A and neutron single particle levels for light to heavy nuclei are compared in figs. |^ and 
Two major differences between the nucleonic and the A spectrum are detected: 

1. A and neutron single particle spectra are overall related by a constant shift and an 
additional quenching factor because the A central potential has a depth of only about 
-30 MeV, compared to -70 MeV for the neutrons (see also fig. ^ and tab. ||). 

2. The spin-orbit splitting of the A states is reduced further being less than what is 
expected from the overall quenching of the potential strength. 

The DDRH calculations reproduce the experimentally suggested very small spin-orbit split- 
ting in A hypernuclei, e.g. p7|j38[] , rather well even without an explicit dynamical suppression 
of spin-orbit interactions as e.g. a A — a; tensor coupling which is used in the QMC model 

i- 
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The reason for the small spin-orbit splitting is understood by considering the evolution 
with increasing mass number. In fig. ^ the spin-orbit splitting of A states for several nuclei 
across the whole mass range is shown. It can be seen clearly that the splitting drops for 
higher masses what is to be expected since it has to go to zero in the nuclear matter limit. 
The splitting also drops in the low mass region - now for the reason that the spin-orbit 
doublets approach the continuum threshold and get compressed before one of them or both 
become unbound. This is a remarkable similarity to the situation found in weakly bound 



neutron-rich exotic nuclei 



For exactly this reason the data point of the ^^Ca Ip shell splitting is, though the 
absolute error is not that large, only of little use to further constrain the relative coupling 
Ri^ (see sec. [Ill C|) . Because the low mass hypernuclei are systematically underbound with 



our standard choice for Rg-^^ the spin-orbit splittings for ^^Ca and ^^Oa are determined with 
a readjusted R^=0.542. [|. 

The explanation of the small spin-orbit splitting is actually found in the overlap of the 
A single particle wave functions and the spin-orbit potential. As seen from binding energies 
and confirmed by the r.m.s. radii of A states shown in tab. |IT1| they are much less localized 
than nucleonic states in the same nucleus. The spin-orbit potential, on the other hand, has 
its strongest contribution always in a peak structure at the nuclear surface, as seen in fig. |^. 
Hence, the overlap of the de-localized A wave functions and the rather sharply localized 
spin-orbit potential is much less than for the stronger bound nucleonic wave functions. As 
a result, a much smaller overlap integral and a reduced spin-orbit interaction energy is 
obtained for A states. The effect of the A spin-orbit potential is largest in the low angular 
momentum doublets of small hypernuclei. Therefore the ideal nuclei to observe spin-orbit 
splitting effects experimentally are those that are heavy enough to bind the Ip doublet just 
without a too large 'close-to-threshold-squeezing'. This can be seen from fig. |^ to be the 
region around calcium. 

D. Systematics of single A states 

In fig. ^ the DDRH single particle spectra are compared to spectroscopic data from 
(tt^^, K~^) reactions. States in intermediate to high mass nuclei are described fairly well 
by the model, while for masses below about ^^SIa deviations to the experimental data of 
up to 2.5 MeV arise. As discussed in section [IV A|, in this mass region the limits of using 



LDA vertices may be approached. We seem to overestimate systematically the strength of 
the repulsive vector interaction in low mass single A hypernuclei. This tendency already 
becomes apparent in going from ^^Va to ^^SIa- In fig. ^ results of a calculation are included 
in which the vector scahng factor was slightly decreased by about 2% to Rtj=0.542. Fig. |l| 
shows that this value is also located in the valley of the distribution. The A separation 
energies in the light mass nuclei are well reproduced now but the agreement in the high 
mass region, however, would deterioriate. The result indicates the sensitivity of the DDRH 



^With this slightly modified R^; = 0.542 the spectra of the low mass hypernuclei are actually 
nicley described; see sec. IVD and fig. ^ 
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calculations on fine details of the interplay of scalar attraction and vector replusion. Clearly, 
under the present conditions a 2% variation is well within the uncertainties of the model 
and, moreover, may be taken to indicate a typical boundary for the validity of LDA in light 
mass nuclei. 

For the heavier nuclei the microscopic DDRH results are of a quality which is at least 
comparable to the phenomenological descriptions. This we consider as a remarkable success 
for a model which essentially contains only a single free parameter, namely the overall vector 
vertex scaling factor R^. 

Finally, we will have a look on [K^j-k^) data measured at the end of the seventies at 
CERN PT| , P2| . This data set was not included in the determination of R^ because of its 
relatively large error bars. Tab. fV\ shows these data together with the DDRH predictions and 
phenomenological RMF |^ calculations. These data contain besides A single particle states 
also informations on the energies of A-particle netron-hole states. These were calculated 
approximately within the Hartree scheme by keeping a neutron hole at the specified place 
during the iteration. A more realistic calculation would require to perform a complete 
RPA calculation which is, however, at present neither feasible nor worthwhile. A further 
complication in modelling these nuclei are the relatively small mass numbers, as discussed 
already above. Taking into acount the fairly low energy resolution of the data and keeping 
in mind the previously discussed subtleties the {K, n) data are also described satisfactorily 
well on a level comparable to phenomenological RMF models. 



V. SUMMARY AND CONCLUSION 

The DDRH theory introduced previously for neutron and proton isospin nuclei was ex- 
tended to hypernuclei by including the full set of SU(3) / octett baryons. Interactions were 
described by a model Lagrangian including strangeness-neutral scalar and vector meson 
fields of qq (q=u,d) and ss quark character. The medium dependence of interactions was 
described by meson-baryon vertices chosen as functionals of the baryon field operators. 
Following DDRH theory their structure is determined such that interaction diagrams con- 
tributing non-perturbatively to the ground state energies and wave functions are cancelled. 
Here, the DDRH vertices were chosen to cancel Dirac-Brueckner ground state correlations. 
Hence, the approach corresponds to a resummation of ladder diagrams into the vertices un- 
der the constraint that infinite matter ground state self-energies and total binding energies 
are reproduced. As the central theoretical result it was found that the structure of Dirac- 
Brueckner interactions strongly indicates that the ratio of nucleon and hyperon in-medium 
vertices should be determined already by the ratio of the corresponding free space coupling 
constants being affected only weakly by the background medium. Apparently, the presently 
available hypernuclear data are, at least, not contradicting such a scaling law. 

Dynamical scaling will have several important consequences for hypernuclear investiga- 
tions. First of all, it might be considered to give a sound theoretical support to the general 
conviction that hypernuclei are suitable to gain information on octett dynamics. A con- 
clusion of equal importance is that hypernuclei should follow essentially the same rules as 
isospin nuclei except for a shift of energy scales. A point worthwhile to be investigated in 
more detail in the future is that hypernuclear scaling might provide a way to study the 
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dynamical evolution of SU(3) flavour symmetry in a medium. 

In order to test dynamical scaling RMF calculations for single A hypernuclei were per- 
formed. Using the previously derived meson-nucleon vertices |T^ and flxing the uA vertex 
by a theoretical value from the literature [0,0] the uo meson-A vertex scaling factor was 
determined by a least square flt procedure thus determining the only free parameter of the 
model empirically from a selected set of data. Calculations over the full range of known sin- 
gle A nuclei led to a very satisfactory description of A separation energies. The deviations 
from the overall agreement for masses below A^16 are probably related to the enhance- 
ment of surface effects in light nuclei. Very likely, they indicate the limitations of static 
RMF calculations with DB vertices obtained in the local density approximation. The minor 
adjustment of parameters, necessary to achieve agreement also for the light mass systems, 
indicates the sensitivity of these surface-dominated nuclei on dynamical details. In a recent 
non-relativistic calculation indeed sizable contributions of hyperon polarization self-energies 
especially in light nuclei were found. 

The results are encouraging and we conclude that DDRH theory, extended to the 
strangeness sector, is in fact an appropriate basis for a microscopic treatment of hypernuclei. 
The present formulation and applications are flrst steps on the way to a more general theory 
of in-medium SU(3)j flavour dynamics. Future progress on dynamical scaling and other 
theoretical aspects of the approach are depending on the availability of Dirac-Brueckner cal- 
culations for the full baryon octett including also the complete pseudoscalar 0~ and vector 
l"*" meson multiplets. 
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TABLES 



TABLE I. Model parameters of the density dependent A-nuclcon model 









i = N 


i = A 


rriN 


939.0 MeV 




26.027 


6.249 


TUA 


1115.0 MeV 


9lii{P=Po) 


6.781 


1.628 


47r 


mo- 


550.0 MeV 


fl2..(p=0) 


40.240 


12.287 


rn^ 


782.6 MeV 


9lii{p=Po) 
Air 


9.899 


3.022 


rup 


770.0 MeV 


9%i 
At: 


1.298 


0.000 



TABLE IL Comparison of the A and neutron central and central rearrangement potential 
depths for single A hypernuclei 





A 


N 


A 




centr. 
c.rearr. 


-3L7 MeV 
1.4 MeV 


-81.5 MeV 
9.3 MeV 


-31.8 MeV 
2.0 MeV 


-85.3 MeV 
8.6 MeV 


S 


-30.3 MeV 


-72.2 MeV 


-29.8 MeV 


-76.7 MeV 






A 


N 


''Oa 

A 




centr. 
c.rearr. 


-31.6 MeV 
2.3 MeV 


-90.6 MeV 
10.1 MeV 


-28.9 MeV 
2.5 MeV 


-86.5 MeV 
8.6 MeV 


S 


-29.3 MeV 


-80.5 MeV 


-26.4 MeV 


-77.9 MeV 



TABLE in. r.m.s. radii of the first orbital momentum states for As, neutrons and protons in 
^OCaA and ^ospbA 





A 


''CaA 

n 


P 


A 


n 


P 


lSl/2 


2.8 fm 


2.3 fm 


2.4 fm 


4.1 fm 


3.8 fm 


3.9 fm 


lP3/2 


3.5 fm 


3.0 fm 


3.0 fm 


4.8 fm 


4.5 fm 


4.6 fm 


lPl/2 


3.6 fm 


3.0 fm 


3.0 fm 


4.7 fm 


4.4 fm 


4.5 fm 


ld5/2 


4.7 fm 


3.5 fm 


3.6 fm 


5.3 fm 


5.0 fm 


5.1 fm 


ld3/2 


6.3 fm 


3.6 fm 


3.7 fm 


5.2 fm 


4.9 fm 


5.0 fm 
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TABLE IV. Transition energies for {K,tt) reactions on a nucleus [ ^Ji^ . These states include 
A particle-n hole excitations of the single A hypernuclei. The experimental values (exp.) are 
compared to DDRH and a phenomenological RMF model [Q] (phen. RMF) with nonlinear a self 
interactions. 









exp 


JJJJrtrl 


pnen.xiivir 




n val. shell 


state 


[ivie V J 


[ivie V J 


[ivie V J 




1^3/2 


(ISi /oA, Ips/olT'^^) 


6.72±2 


6.69 


5.02 






(IP3/2A, lp3/2n~^) 


18.48±2 


15.11 


17.21 


''Oa 


lPl/2 


(IS1/2A, lpi/2n-'^) 


3.35±2 


5.76 


3.53 






(IS1/2A, lp3/2n"^) 


9.90±2 


10.13 


9.46 






(IP1/2A, lpi/2n-i) 


13.20±2 


16.16 


13.89 






(IP3/2A, l?'3/2'^~^) 


19.20±2 


18.40 


18.88 


''CaA 


14/2 


(IP1/2A, ld3/2n-i) 


5.79±2 


8.84 


7.40 






(I4/2A, ld3/2n-i) 


14.47±2 


11.34 


15.48 






(W5/2A, l(i5/2n"^) 


19.35±2 


20.07 


20.71 
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FIGURES 



FIG. 1. distribution for the deviation of DDRH A single particle energies and hypernuclear 



data (obtained in (tt, K) reactions [32-3g]). In the calculations, the scalar and vector vertex factors 
[Ra,Rw) were varied freely. 
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FIG. 2. The radial variation of the DDRH A meson vertices Ta,uj (upper graph) and of the 
vector densities (lower graph) for lsx/2 ^ states in light and heavy nuclei. Results for the single A 
hypernuclei ^'^^PbA (full), ^^Oa (short dashed) and ^Bca (long dashed) are displayed. 
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FIG. 3. Comparison of the Schrocdinger-equivalent A and neutron central potentials including 
rearrangement. Results for ^'^'^PbA and ^^Oa are shown in the upper and lower part of the figure, 
respectively. The A and neutron potentials are of a similar shape but different depth. 
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FIG. 4. Mean-field 

potentials for single A hypernuclei. The lowest order Sehroedinger-equivalent conventional po- 
tential (solid line) and rearrangement potential (dashed line) are shown. It is clearly seen that the 
polarization effects described effectively by the rearrangement potential are most important for 
light nuclei. 
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FIG. 5. DDRH A single particle spectra in light to heavy hypernuclei. 
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FIG. 6. DDRH neutron single particle spectra in light to heavy hypernuclei. 
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FIG. 7. Dependence of the spin-orbit splitting for lp3/2,i/2, ld5/2,3/2j lf7/2,5/2, lg9/2,7/2 and 



lhxi/2,9/2 ^ states on nuclear mass. The only available datapoint |37] for the lp3/2,i/2 doublet in 
-"^^Ca is also shown (see sect. QIC). 
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FIG. 8. A and neutron spin-orbit potentials in ^'^^PbA and CaA. The lowest order 

Schrocdinger-cqnivalont potentials are shown. 
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FIG. 9. Comparison of DDRH and experimental single particle energies. The lines are drawn 
to guide the eye. Results for R^=0.542 are indicated by a dashed line. In order to remove finite 
size and, especially, surface effects, the energies are shown as a function of A~^/'^. For A— > they 
converge asymptotically to the binding energy of a single A in infinite matter, Ea=-28 MeV. The 



data originate from X(tt, K) X\ reactions [32-36| 
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